Optical conductivity from local anharmonic phonons 
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Recently there has been paid much attention to phenomena caused by local anharmonic vibrations 
of the guest ions encapsulated in polyhedral cages of materials such as pyrochlore oxides, filled 
skutterdites and clathrates. We theoretically investigate the optical conductivity solely due to these 
so-called rattling phonons in a one-dimensional anharmonic potential model. The dipole interaction 
of the guest ions with electric fields induces excitations expressed as transitions among vibrational 
states with non-equally spaced energies, resulting in a natural line broadening and a shift of the 
peak frequency as anharmonic efi'ects. In the case of a single well potential, a softening of the 
peak frequency and an asymmetric narrowing of the line width with decreasing temperature are 
understood as a shift of the spectral weight to lower level transitions. On the other hand, the 
case of a double minima potential leads to a multi-splitting of a spectral peak in the conductivity 
spectrum with decreasing temperature. 

PACS numbers: 63.20.Ry, 63.20.Pw, 78.20.Bh 
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I. INTRODUCTION 

Anharmonicity in lattice vibrations has been one of 
the old problems in condensed-matter physics^'^. An- 
harmonic effects in acoustic phonons were treated by per- 
turbation theory^i^, while those in local vibrations were 
investigated in impurities or disordered systemsii^i^. It 
was pointed out that the effects of the anharmonicity 
in local vibrations appear in the characteristic tempera- 
ture dependence of the vibrational frequency and of the 
line width^. Since an isolated irregular atom receives an 
anharmonic potential from the regular lattice, analysis 
have been made mostly on the anharmonic oscillation re- 
ceiving effects of surrounding oscillation of the regular 
lattice. 

Recently a revised interest on anharmonic phonons 
has arisen in relation to a material series of 
pyrochlore oxides§i§,^ filled skutterditesi^iiiii^ and 
clathratesi^ii^ii^'i^iii. Those materials, which are usu- 
ally electrical conductors like metals, semimetals or 
heavily-doped semiconductors, have a common feature 
that some numbers of atoms form a three-dimensional 
network of polyhedral cages, in each of which a guest 
ions is accommodated. When the cages are oversized, 
the guest ion vibrates with a large amplitude in an an- 
harmonic potential. Such vibrations are named as the 
rattling phonons. To note, depending on kinds of guest 
ions, on-centering or off-centering vibrations occur even 
in the same cage structure. 

There have been reported various anomalous phenom- 
ena in the above materials, some of which have been dis- 
cussed in relation to those rattling phonons^^. In apply- 
ing some clathrate compounds to thermoelectric material 
devices, for example, the rattling phonons, in particu- 



lar off-centered, are expected to suppress strongly the 
thermal conductivity by effectively scattering Debye-like 
acoustic phonons propagating through the cage network 
and carrying heat entrop y^^i^^'^^i^° . An alternative ex- 
ample is found in the superconductivity in a /3-pyrochlore 
oxide KOS2O6. It was suggested that rattling vibrations 
of the K+ in the OsOg octahedral cage were responsi- 
ble for the strong-coupling superconductivity and also 
for an electron-mass enhancemen t^^i^^ . It may be fair, 
however, to state that these interesting issues as for the 
question how rattling phonons interact with cage acous- 
tic phonons and/or charge carriers are far from being well 
understoodiSi. 

So far lattice vibrational modes including rattling 
phonons in the above cage materials have been stud- 
ied rather extensively with use of spectroscopic measure- 
ments such as inelastic neutron and Raman scatterings^^. 
Some low-lying rattling modes are clarified to exhibit 
softening with decreasing temperatures. The softening 
phenomenon has been well recognized as one of the an- 
harmonic effects from rattling phonons, which was dis- 
cussed with a quasi-harmonic approximation^''. Beside 
these spectroscopies, an infrared-active optical measure- 
ments, particularly in the Terahertz range, would be, in 
principle, a powerful tool to clarify the charge dynamics 
in low-lying optical phonons near g ~ with available 
optical conductivity spectra. 

Recently, time-domain terahertz spectroscopy^^ has 
been successfully applied for the first time to observe 
the rattling phonons around ITHz in a type-I clathrate 
BasGaieGcso (BGG)^. In this paper, we investigate 
systematically the optical conductivity spectra from the 
rattling phonons in an on-centered and off-centered po- 
tential, and apply the obtained theoretical classifica- 
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tion of the spectra to analyse the experimental result 
of BGG. The detail of the experiment will be presented 
elsewhere^^. We express the optical conductivity by level 
transitions among states in an anharmonic potential^. 
The present analysis shows that the most important ef- 
fect of an anharmonic potential is the non-equal energy 
spacing of level transitions, resulting in an intrinsic and 
asymmetric spread of the excitation spectral, in contrast 
to the harmonic case of the equal energy spacing simply 
resulting in the Lorentzian spectral shape. Also there 
arise variety of level schemes for low-lying states depend- 
ing on an on-centered or off-centered potential. Those 
features induce the natural broadening of the line width 
at higher temperature, softening of the peak frequency^ 
and multi-peak structures at low temperature in optical 
conductivity. 

In the present paper, a one-dimensional model is used 
for simplicity. Essential effects from the anharmonicity 
are included. In the next section, the model and the 
expression of the optical conductivity are presented. In 
Sect. 3 some of numerical results are presented. Compar- 
ison with the recent experiment22. is discussed. Sect. 4 is 
devoted to the conclusion. 



II. FORMULATION 

As a model to describe the motion of a guest ion in the 
cage, we take the following one-dimensional anharmonic 
potential model, for simplicity. 
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(2.1) 



where M, p and x are the mass, momentum and spa- 
tial coordinate of the guest ion, respectively. We neglect 
effects of acoustic phonons and electrons. 

The optical conductivity from the guest ion is obtained 
by considering the polarization induced by an applied 
oscillatory electric field E{t), 



Hi = -qE{t)a 



(2.2) 



where q is the charge of the guest ion. By use of the 
linear response theory, the polarization P{t){= {qx{t))) 
is obtained as 

P(t) = ^J dt'q^{Rx(t)x{t'))E(t') , (2.3) 

where (• • • ) indicates the thermal average and "i?" means 
the retarded function. By taking the Fourier transform, 
the polarizability a{u}) defined by 



P{uj) = a{uj)E{uj) 



is obtained as 



a{uj) 



q^N 



(2.4) 



(2.5) 



where the density of the quest ion N is taken into account 
and Gxx('^) is defined by 

{Rx{t)x{t')) = ^ y dc^e-^"(*-*')G,,(a;) . (2.6) 

Let us denote eigenstates and eigenvalues of the Hamil- 
tonian i/ by |rt) and i?„, respectively. 



H\n) ~ En\n) 
Then Gxx{^) is expressed as 



(2.7) 



Gxx{i^) = ^ |(n|a;|m)|^ 



)/Z 



io + iTo/2-iEn- Em)/h 

(2.8) 

with P = 1/kBT and Z = e . Here we have intro- 
duced a phenomenological parameter of the decay width 
Fq. Since the optical conductivity and the polarizability 
are related to each other as 

a{uj) — ~iuja{uj) , (2-9) 

we have the complex optical conductivity as 

a{uj)/(TQ 



^ \{n\{x/x,)\m)\' 



1 



with 



LO - UJn7n+i'^o/2 + CJ„m + /2 



{En - E^)/h . 



(2.10) 



(2.11) 



Hereafter, we use the notation E„rm Wnm and i^nm, re- 
spectively, as the energy, angular frequency and fre- 
quency for the transition from the m-state to the n-state. 
The normalization for the conductivity ctq is given by 



ao = q^Nxl/n . 



(2.12) 



with a length scale xq determined shortly. In the follow- 
ing analysis, we choose a suitable energy scale Tiwo, which 
can be, for example, the energy corresponding to ITHz 
or a harmonic frequency evaluated from the quadratic 
term in the potential, or an observed phonon energy. The 
scaled parameters are defined as 



and the length scale xq is given as 



Xo 



Mujf) 



(2.13) 



(2.14) 



In the next section, we will calculate <j{uj) for cases 
of fc > 0, fc = and A: < 0, and discuss characteris- 
tic features of the optical conductivity from the rattling 
phonon. 
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III. NUMERICAL ANALYSIS 

A. The case of fc > 

We can choose k = 1 witho ut los s of generahty, In 
this case, ujq is given by ujq = y^k/M, the harmonic fre- 
quency. Let us estimate the magnitude of A. In the 
band structure calculation for BagGaigGescr^, the coef- 
ficients k and A were evaluated as ka\/2 = 2.813mRy 
and Aa^/4 — 1.602mRy, where is the Bohr ra- 
dius. Then using the mass of the Ba-atom 137.35 a.u., 
we have the harmonic frequency vq — 0.6975THz and 
A = 4.292 X 10~^. The actually observed phonon fre- 
quency lowest lying in BasGaieGeao is about ITHz and 
the parameters k may vary about twice or so. Taking into 
account the above estimation and Eq. (|2.13p . we take the 
value of A in the range of A = 0.0 — 0.05. 
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FIG. 1: (Color online) (a) Optical conductivity for various A 
with = 1 at kBT/fiuJo — 10.0, (b) excitation energy ujio and 
(c)energy eigenvalues for A = 0.05. 

In Fig. [1] (a) we plot the optical conductivity for var- 
ious A at kBT/fiujQ — 10.0. Hereafter, we introduce a 
constant decay width ro/2u;o = 0.1 by hand, in order to 
smooth the frequency dependence of the optical conduc- 
tivity. How each level transition has a decay width de- 
pends on interactions with electrons or acoustic phonons. 
Details of such effects are not considered in this paper. 
Also in the following, denotes the real part of the 

complex conductivity. At A = 0, the phonon mode has 
the single angular frequency oj — uj{), and the line width 
is the decay width Fq. As A increases, both of the peak 
position and the line width increase. That is, the line 
width at higher temperature is determined by the anhar- 
monic parameter A. In Fig. [T] (b), we plot the excita- 
tion energy W(„_|_i) „ = — En)/h. This is obtained 
by calculating the energy eigenvalues numerically. One 



example of the behavior of energy eigenvalues is shown 
in Fig. [He) for A = 0.05. The label "n" is identical 
with the boson number in the harmonic case of A = 0. 
The energy eigenvalues deviate from the linear behavior 
as En ~ Co -l- cin + e2fi^ but a perturbation calcula- 
tion is not applicable, since Xn? becomes an order of one 
for n ^ 10. The calculation shows that the transition 
probability arises mostly from (n -t- l|a;|n). As is seen in 
Fig. [TJb), the excitation energies increase as eigenval- 
ues(i.e. n) increase due to the anharmonicity; the larger 
the A, the more spreading excitation energies become. 
The non-equal energy spacing of the phononic level tran- 
sitions leads to the intrinsic spread of the line width in 
the optical conductivity, as is seen in Fig. [Ija). It should 
be noted that excitations with larger n ('~ 30) contribute 
at kBT/fwjQ = 10. 
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FIG. 2: (Color on line)(a)Temperature dependence of the op- 
tical conductivity for A = 0.04 and k = 1. (b) The corre- 
sponding spectral weight at fc_BT/?iaJo=0.05,5.0 and 10.0. 

In Fig. [IJa), we plot the temperature dependence of 
the optical conductivity. The optical conductivity shows 
the softening of the peak position and the narrowing of 
the line width with decreasing temperature. At low tem- 
perature, only the first level transition contributes. As 
temperature increases, the transition involves higher en- 
ergy levels, resulting in a shift of the effective peak posi- 
tion. Also the line width increases, since more level tran- 
sitions contribute. In order to see the line-distribution 
without the smoothening by the width Fq, we plot, in 
Fig. IHb), the spectral weight in the optical conductivity 
given by 



F{uj)/ao = TT 



E 



|(n|e|m)p 



-I3E„ 



(3.1) 



for kgT/llujQ = 0.05,5.0 and 10.0. Sharp lines arising 
from the non-equal spacing of the level transitions dis- 
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tribute rather densely in the present parameter values. 
Because of the Boltzmann factor and a factor n from 
|(n|^|n ± the maximum position of the weight in- 
creases with temperature and decreases exponentially as 
the energy hu; increases. At T — OK, the main con- 
tribution arises from the transition between the ground 
state and the first excited state. The softening of the 
peak frequency and the sharpening of the line width with 
decreasing temperature are direct consequences derived 
solely from an inequivalence of the level spacing that is 
essential in the anharmonicity. Therefore measurements 
on the temperature dependence of the optical conductiv- 
ity can provide a direct evidence for the anharmonicity. 
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FIG. 3: A-dependence of uio for k = 1 

In Fig. [3] the transition energy for the 0-^1 transition 
is shown, which gives the peak frequency at T = OK. In 
the present parameter region, the curve is fitted by 

wio/wo = 1 + 0.7414A - 0.7998P , (3.2) 

which deviates above A ^ 0.05 from the result of the 
perturbation, wio = 1 + (3/4)A - {9/8)X'^ (dashed line). 



by red solid lines, upper and lower energies which give 
half values of the peak intensities. The red dashed line 
is for their averaged energy (say, mid- frequency) . The 
blue diagonal-crosses indicate the result obtained from 
the following equation in the mean field theory27 (The 
factor 3 is corrected.). 



1 + 3X 



1 



exp 



hcjQ ijj 



+ 



- 1 



\ 



J 



(3.3) 



We see that the averaged value (red dashed line) well 
agree with the mean field results (blue diagonal-crosses). 
However, our result of the peak structure is asymmetric, 
reflecting the spectral weight, Eq. p.ip . We express the 
line width at temperature T as 



r(r) = 7„(T) + 7,(T) 



(3.4) 



with 7m (r) (71 (r)) being the upper (lower) half width 
from the peak frequency, which is shown in Fig. |4] (b) . 
The saturation of the line width is due to the decay width 
Fq, which may originate from a phonon-electron interac- 
tion and a rattling phonon-acoustic phonon interaction 
neglected in this paper. However, at higher temperature, 
r(r) increases fairly larger than Fq. Also the broaden- 
ing in upper and lower frequency is asymmetric. The 
line broadening due to the anharmonicity leads to a non- 
Lorentzian spectral shape. 

From the analysis of this subsection, we see that an- 
harmonicity of the rattling phonon leads to softening of 
the phonon frequency and sharpening and asymmetric 
change of the line width with decreasing temperature. 
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FIG. 4: (Color online) Temperature dependence of (a) the 
peak position and (b) line width for A = 0.04 and k = 1 



FIG. 5: (a) Energy eigenvalues and (b) excitation energy for 
A: = 



In Fig. 3] (a) we plot the temperature dependence of 
the peak frequency for Fo/2u;o = 0.1 (red circles) and. 



When fc = (or |A:| << Aa^), the peak position and 
the line width are solely determined by A. We choose 
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huo as a characteristic energy of the system, for example, 
a;o/27r=lTHz. In Figs. \E\{a) and (b), the energy eigenval- 
ues and the excitation energy W(„+i).„ = — En)/h 
are plotted. The lower excitation energies are softened 
and, just as well as in the case of fc > 0, the effectively 
contributing n decreases with decreasing temperatures, 
resulting in the characteristic behavior, softening and 
sharpening, of the optical conductivity as is shown in 
Figs. [6] and [T] Note that the lowest excitation energy is 
softened up to about a half of the characteristic energy 
because of the shallow shape of the potential bottom. 
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FIG. 6; (Color online) Temperature dependence of the optical 
conductivity for A = 0.04 and k — 0. 
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FIG. 7: (Color online) Temperature dependence of the optical 
conductivity for A = 0.1 and k — 0. 

In the case of A = 0.04 (Fig. [6]), the difference of 
the neighboring excitation energies is less than Fo, so 
that the optical conductivity shows a smooth softening, 
though the shift is enlarged and some shoulder is seen 
at low temperature (fcsT/Tiwo = 0.5). Such a shoulder 
structure is smoothed for a smaller A. When A becomes 
large, the difference of the lower excitation energies ex- 
ceeds Fq, and as in Fig. [7] the optical conductivity shows 



a double peak structure as temperature decreased. The 
upper peak is reduced as temperature is further lowered. 
The broadening of the line width at higher temperature 
is enhanced as A increases. 



C. The case of fc < 

For negative fc, the anharmonic potential has double 
minima. We investigate the temperature dependence of 
the optical conductivity for A — 0.04 by changing k. De- 
pending on the depth of the double minima, various pat- 
terns of temperature-dependent optical conductivity are 
obtained. We choose hojQ as a characteristic energy of 
the system, for example, tJo/27r=lTHz. 

In Fig. [51 we show the temperature dependence of the 
optical conductivity for various fc. The softening of the 
peak position is seen in Figs. [8Ua)-(c), but in Figs.[8l^d)- 
(f) the softening of the peak frequency is followed by the 
hardening, with decreasing temperature. 

Low temperature behaviors have various variety, es- 
pecially there appear structures with double or triple 
peaks, though it depends on the magnitude of the width 
Fq. This can be understood from the low energy level 
transition in the double well potential. In Fig. [H^a) , one 
example of the potential with fc < is illustrated to- 
gether with the eigenenergy levels. Low energy levels are 
much modified by the depth of the potential well. In 
Fig. [5Jb) the fc-dependence of En for low energy lev- 
els are plotted. Dotted lines are the potential double 
minima with the value V {±^y—k/X) / hwo = — fc^/4A and 
the local maximum V{Q) = 0, respectively. When — fc 
increases, the levels {Eo,Ei), {E2,Ej,), successively 
degenerate forming low- lying tunneling modes. In Fig. 
|9l[c) the transition energies w„m are plotted among low- 
lying eigenstates. The dotted line is for the depth of the 
potential well. Transition energies wio and ^32 show suc- 
cessively softening as — fc increases, while the neighboring 
excitations uj2i^ i^iz increase and become larger than ex- 
citation energies among higher levels (o.'ii.io in Fig. EI^c)). 
The upturn of the excitation energies of LO21 is correlated 
with the depth of the potential well. In the following we 
discuss more details of low temperature behaviors in Fig. 

m 

Figs. El^a) and (b) are for shallow double wells. Tran- 
sition energy of huiQ is softened and it becomes smaller 
than 111021 due to the effect of the double well. Further, 
{lj2i — ^w) is larger than Fq, so that there appears two 
peaks at low temperature and the higher peak of LO21 di- 
minishes as T — > OK. The peak around lo/ujo = 1 which 
remains even at T = OK corresponds to the transition 

Fig. [8ljc) is for fc = —0.3. The states and 1 become 
very close but still have finite difference. We can identify 
lower peak as wio, LO21 and W30. The peak for W21 is 
reduced with decreasing temperature (See also Fig. [H^c)). 

Fig. [Hid) is for fc = -0.4. The states and 1 are 
almost degenerate, and this soft mode does not appear 
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FIG. 8: (Color online, Wide figure) Temperature dependence of the optical conductivity for various A; < 0. Parameters 
areA = 0.04, ro/2cJo ~ 0.1. t is the reduced temperature t — kBT/tiLJo. 




FIG. 9: (Color online) (a)Schematics of the double-minima 
potential and energy levels for k = —0.4, and fc-dependence 
of {h)En and (c)aj„mfor low energy levels. A is chosen as 0.04. 



in the optical conductivity because of the factor lu in Eq. 
(|2.10p . The state 2 is inside the double weU and the state 
3 is above the maximum at a; = 0. Then 0^21 is larger 
than 0^32. We can identify the peaks in Fig. [SJd) as 
^32, ^12 and UJ03 from the low energy side. Note that 
the intensity of the peak UJ32 increases and decreases as 
temperature decreases. 

In Fig. ^e), k is further reduced as A; = —0.5. The 
state 2 and 3 start to degenerate and W32 becomes lovifer 
than a;2i- We can identify the peaks in Fig. [8{e) from 
the low energy side as 0^43 and (a;i2,iu;o3). 

In Fig. ^f), the state 2 and 3 are completely degen- 



erate and cl;2i is the main transition at low temperature. 
Since this energy is much larger than energies for higher 
level transition, the peak position first decrease at higher 
temperature and then increase at low temperature. 

When k is further reduced and the double minima be- 
come deep enough, the hardening of the peak frequency is 
obtained rather than softening as was discussed in ref.— . 

In this way, we have various patterns depending on the 
strength of the double well 



D. Comparison with the experiment 

We have performed the time-domain THz spectroscopy 
and obtained the optical conductivity^^ in a type-I 
clathrate B&sGaieGesir^'^^^^^. The details of the 
experimental method and analysis will be presented 
in a separated paper—. The obtained temperature- 
dependence of the phonon spectral for the lowest mode(~ 
1.2TIIz) is shown by colored symbols in Fig. (TUl No 
multi-peak structure is observed. Then from the pat- 
terns presented in this paper, we conclude fc > and we 
take k = 1. 

We adjust coq, A, Tq and cto to fit overall behaviors and 
specially the higher frequency region (1.25^1.40 THz) 
of the spectral line. We choose uj^jl-n — 1.143THz, A — 
9.66 X 10-3, ro/27r = 0.067THz, and ctq = % .%2Qr^ cix\~^ . 
Theoretical results are presented by solid lines in Fig. [TOl 
The agreement between the experimental and theoretical 
results in temperature-dependence is very good. 
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FIG. 10: (Color online) Comparison of the optical conductiv- 
ity for BagGaieGeso 
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FIG. 11: (Color online) Temperature-dependence of the peak 
frequency and the line width 



The temperature dependence of the peak frequency 
and the line width are estimated from the experimental 
data by the Lorentzian fit and are plotted by the solid 
rectangulars in Fig. [TTTa) and (b), respectively. Theoret- 
ical results are plotted by solid and dashed lines. In Fig. 
Illf a). the red solid line is for the peak frequency I'poak 
and the green dashed line is for the mid- frequency I'mid- 
In Fig. [TIT b). the red solid, blue dashed and light blue 
short-dashed lines are frequencies ^'r/2, v^^ and v^^ cor- 
responding to the half width, the lower and upper half 
width, respectively. The experimental peak frequency 
situates between the theoretical peak and mid-frequency. 
Also the half-width is between the theoretical half-width 
and lower half-width. Since the Lorentzian fit is apt to 



lead a higher peak frequency and a narrower line width 
in an asymmetric line shape, the agreement between the 
experimental and theoretical results is reasonably good. 
The present comparison shows that the temperature- 
dependence of the rattling phonon in BagGaigGeao is well 
described by the anharmonicity of the potential without 
considering details of effects from other interactions. 

It has been reported that a type-I clathrate 
BasGaieSnao has an off-centered potential of the guest 
ionii. It is an interesting problem to see if the optical 
conductivity of this material shows a behavior with fc < 
discussed in this paper—. 



IV. CONCLUSION 

In this paper we have investigated theoretically the 
temperature dependence of the optical conductivity from 
the rattling phonon. The guest ion feels an anharmonic 
potential from the cage, and the anharmonic effect ap- 
pears characteristically in the temperature dependence, 
that is, the softening of the peak frequency and sharpen- 
ing of the line width with decreasing temperature. 

In the case of the positive quadratic term, one can ex- 
pect the quadratic coefficient is roughly determined from 
the saturated peak frequency at low temperature and the 
quartic term is determined from the line width and shift 
of the peak frequency. 

In the case of the negative quadratic term, various pat- 
terns of the optical conductivity are expected depending 
on the strength of the double minima in the potential. 
Multi-peak structures appear at low temperature and in- 
crease and decrease of the peak frequency with temper- 
ature are obtained. 

We have shown that measurements on the tempera- 
ture dependence of the optical conductivity can provide 
a direct evidence for the anharmonicity. 
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